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1. Introduction 

The Nakajima-Zwanzig projection method [1, 2] is central to modern open 
quantum systems theory [3, 4]. The method provides a systematic theoretical 
approach to Markovian dynamics and sets a basis for a systematic (perturba- 
tive) approach to non-Markovian dynamics [3, 4]. The open systems theory 
provides the foundations of quantum measurement theory [6, 7], decoherence 
[6, 7, 8], and the emergence of thermodynamic behavior [9]. Applications of 
the theory of open quantum systems are found in practically all areas of 
physics, ranging from quantum optics [10], quantum information [11] and 
condensed-matter physics [12] to chemical physics [13] and spintronics [14]. 

The key idea behind the Nakajima-Zwanzig projection method consists 
of the introduction of a certain projection operator, V, which acts on the 
operators of the state space of the total system "system+environment" (S + 
E). If p is the density matrix of the total system, the projection Vp (the 
"relevant part" of the total density matrix) serves to represent a simplified 
effective description through a reduced set of variables. The complementary 
part (the "irrelevant part" of the total density matrix), Qp — (I — V)p. 
For the "relevant part", Vp(t), one derives closed equations of motion in 
the form of integro-differential equation. The open system's density matrix 
ps(t) = tr E Vp(t) contains all necessary information about the "relevant" 
degrees of freedom of the open system S. 

In this paper we point out a caveat regarding applicability of the Nakajima- 
Zwanzig projection method. This subtle observation sets a limit on appli- 
cability of the method for certain tasks in the open systems theory. This 
observation is the true subject of the present paper. Here we just want to 
give a clue, which can be expressed through two mutually related questions: 
Is it allowed to project the total state in regard to the different sets of the 
degrees of freedom, i.e. to the different system-environment splits? If it is 
allowed, can it be simultaneously performed? The answers provided below 
read, respectively, "yes" and "no". 

However, the universally valid quantum mechanics does not support the 
later. There is not any a priori privileged system-environment split [15-21]. 
Quantum decoherence is typically studied starting from a fairly unprincipled 
choice of system-environment split. In this sense, decoherence is often con- 
sidered as an approximate process [7, 8, 22]. In practice, division of the total 
system into subsystems (into "(open) system" and "environment") is "im- 
perfect" [7] but potentially useful-e.g. one can try to devise the more useful 
recipes for drawing quantum correlations from a composite system [17, 20]. 

So we arrive at the main result of this paper: the method of projecting 
state of the total system typically fails in simultaneously describing the alter- 
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native decompositions (structures) of the total system "system+environment" . 
In other words, the Nakajima-Zwanzig projection method cannot be used to 
simultaneously describe dynamics of the different decompositions (partitions 
in to subsystems) of the total system. 

2. Simultaneous dynamics of the structures 

We are interested in a composite system C that is described by a bipartite 
structure C = S + E. Providing the open system S is described by a proper 
master equation, we wonder if the knowledge of the S"s dynamics can help in 
describing dynamics of an alternative open system S', which defines an alter- 
native structure of the composite system C; C = S' + E'. While description 
of the different structures at the same time is basic-notice the simultaneous 
redefinition of both the "system" and "environment" -it can also have some 
practical motivations. E.g. we can wonder if the S' system can be more eas- 
ily accessible in a laboratory than the S system. This can provide the more 
convenient recipes for manipulating the C's degrees of freedom. Or we may 
be interested in dynamics of the S' system, which is not directly accessible 
in a laboratory. 

Quantum mechanics is insensitive to different structures (decompositions 
into parts) of a composite system C. That is, quantum mechanics equally 
treats the different structures of C. The von Neumann- Liouville equation 
(R = 1): 

^ = -[**)] ID 

where H is the total system's Hamiltonian acting on the total system's 
Hilbert state-space H, equally applies to every decomposition (structure) 
ofC. 

Consider a pair of structures, S + E and S' + E'; S + E = C = S' + E'. 
This provides the different tensor-product-structures for the total Hilbert 
state space, Tig ® We = H = Us' ® He 1 - Of course, the total system's 
Hamiltonian and state, as well as reduced state of any subsystem (obtained 
by the proper tracing out) are unique in every instant of time. The different 
forms of the total Hamiltonian 



H (se) = H S + H E + H SE = H = H s , + H E , + H s >e> = H^ S ' E '\ (2) 

where the double subscripts distinguish the interaction terms. Then eqs. (1) 
and (2) provide simultaneous description of the reduced dynamics for both 
open systems, S and S': 
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^ = -it ri [H,p(t)}, t = S,S', (3) 

where pi(t) = trip(t); for the i = S, the Hamiltonian H takes the form H^ SE \ 
while for the i = S', the Hamiltonian takes the form H^ s E \ 

The technical difficulties in solving equations eq. (3) have historically 
led to the development of the different methods, notably to the Nakajima- 
Zwanzig projection method [1-4], which introduces a projection operator V 
and its complementary projection operator, Q = I—V. The projection Vp(t) 
is required to contain all necessary information about the open system S: 



Ps(t)=tr E p(t)=tr E Vp(t),Vt. (4) 

Then the task is to provide a closed master equation for ps(t), such as 
e.g. the generalized Nakajima-Zwanzig equation or the time-convolutionless 
master equation [3, 4]. 

The projections can be defined [3, 4, 8, 23]: (i) Vp(t) = (tr E p(t)) Cg> 
p E [with arbitrary p E ], (ii) Vp(t) = J2n( tr EPsnP(t)) <g> p En [with arbitrary 
orthogonal supports for p E s], and (hi) Vp{t) = J2i(t r EPEiP(t)) ® Pe% [with 
arbitrary projectors for the E system]; by P, we denote the projectors on 
the respective Hilbert state (factor) spaces. 

To illustrate, consider the standard quantum teleportation setup [13] for 
three qubits, C = 1 + 2 + 3. The two bipartite structures, S + E = 1 + (2 + 3) 
and S' + E' = (l + 2) + 3. Then the quantum state of C [13] can be 
written as 



4 1 

\u) s ®\4>)e = I*) =Eo^ s,<8) ^ E '- ( 5 ) 
1=1 1 

The projection (i) gives V'p = (1/4) X)» |*)s'(*| ® Pe> for the S' + E' 
structure that is a mixed state fulfilling the condition Eq. (4). Application 
of the same procedure for the S+E structure, i.e. the projection (i), gives the 
equality V p = p for the S+E structure also satisfying eq. (4). The projector 
V' provides mixed state for the total system, and consequently mixed states 
for both 1 and 2 subsystems. However, the V projector provides a pure 
state for the 1 subsystem-see the lhs of Eq. (5). So, the two projection 
operators, V and V, cannot equal to each other. Furthermore, as they 
provide different states e.g. for the 1 subsystem in an instant in time, the two 
projectors exclude each other. The same conclusion applies to the hydrogen 
atom differently structured either as "electron+proton (e + p)" or as the 
"center of mass + relative position (CM + R)" . The atom's (instantaneous) 
state |$) [19, 24] 
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J2 c i\ l )e ® \l) P = 1$) = \x)cm ® \nlmim s ) R \ (6) 
i 

in Eq. (6), the quantum numbers n,l,mi,m s are the standard quantum 
numbers known from the quantum theory of the hydrogen atom. 

The equalities eqs. (5) and (6) are instances of a quantum mechanical 
rule of "entanglement relativity" (ER) [15, 18-20, 25, 26] (and the references 
therein), which has recently been extended to relativity of the more gen- 
eral non-classical (quantum) correlations quantified by "discord" -quantum 
correlations relativity (QCR) [26]. Both ER and QCR emphasize: transfor- 
mations of the degrees-of-freedom typically effect in a change in correlations 
present in the composite system C . Typically, for arbitrary instant in time, 
a factorized state for one structure is endowed by non-classical correlations 
for the alternate structure, e.g. for a mixed state p (p 2 ^ p, trp = 1) 

ps ® Pe = P = ^Ps'i ® Pm, = L ( 7 ) 

i i 

If eq. (7) is a consequence of the projection (i) (when the p in eq.(7) 
should be exchanged by Vp) then it contains all necessary information about 
the open system S, i.e. the requirement eq. (4) is fulfiled. However, the rhs 
of eq. (7) is in general not of any type of the above distinguished projections 
(i)-(iii). Likewise for eqs. (5) and (6), the projection J2i Kpsh® Pe'% does not 
in general encapsulate all necessary information about the open system S' . 
Rather, the "irrelevant part", Qp, can be expected to bring some information 
about the open system S'. 

Given eq. (4) is fulfilled, eq. (7) (with the Vp instead of p) implies: 

tr E Qp(t) = tr E (p(t) - Vp(t)) = tr E (p(t) - p s (t) ® p E ) = 0, V*. (8) 
The analogues condition regarding the S' + E' structure: 

tr E ,Qp(t) = tr E >(p(t) - p s (t) ® p E ) = 0,Vt. (9) 

Then both tr E Vp(t) = ps(t) and tr E >Vp(t) = ps'(t) and one can write 
sumultaneous master equations for the S and S' systems, with the constraints 
coming from eq. (7). If eq. (7) refers to the initial state, p(t ), then both 
eq.(8) and eq. (9) are fulfilled in the initial instant of time t a , see also eq. 
(32) below. 

However, eq. (9) is typically not fulfilled. More precisely: 
Lemma 1. For the most part of the composite system's dynamics, validity 
of eq. (8) implies nonvalidity of eq. (9), and vice versa. 
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Proof: We refer to the projection (i) in an instant of time: 



Vp = (tr E p) ®p E - (10) 

Similar arguments can be used for the cases (ii) and (iii). 
A) Pure state p = |^)(^|, while, due to eq. (10), tr E Q\^)(^\ = 0. 

We consider the pure state presented in its (not necessarily unique) Schmidt 
form 

|*> = £q|0s|^, (11) 

i 

where ps = trE\^)($\ = HiPi\i)s(i\-> Pi = l c i| 2 and for arbitrary p E ^ 
tr s \^)(^)- Given p E = ir a \a)E((x\, we decompose \^/) as: 

|^} = E c ^>sl«>^ ( 12 ) 

with the constraints: 

j2h\ 2 = i = J2 t , E ic-l 2 = i, vi, (13) 

i a. a 

Then 

= mm - Y,Pi*o\i)s(i\ ® W)e{o\. (14) 

We use ER: 

\i)s\a) E = J2 D mn\™)s'\n) E ' (15) 

m,n 

with the constraints: 

J2 D ^Dit = ^5 aa , (16) 

m,n 

With the use of eqs. (12) and (15), eq. (14) reads: 



E [ E ^W'o^^mV-Eft^^^ 

(17) 



m,m'n,n' i,v ,a,a i,a 



After tracing out, tr^: 



QCjaC*, C*, Q , D^ n D^'* - p^ a D^ n DZr n \ \m) s >(m'\ (18) 
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Hence 



*r^e|*>(*| =0«EE c i C Ut c$C$jD^D™;- mtt D% n D%' n l = 0, Vm,m'. 

i,a,n i',a' 

(19) 

Introducing notation, A™ = Y^i, a c iCi a Dl^ n , one obtains: 
fr^ei^X*! = & A mm , = Y\KK'* -EPi^D^D^] = 0, Vm,m'. 

n i, a 

(20) 

Notice: 

£4»m = 0. (21) 

m 

which is equivalent to trQ|^)(^| = 0, see eq. (14). 
B) Mixed (nonentangled) state. 

P = ^2^iPSiPEi, PSi = y^Pim\Xim)s(Xim\, 
i m 

In eq. (22), having in mind eq. (10), tr E Qp 

and p E = Eg^lV'gMV'.zl ^ tr sP- 
Constraints: 

£ K = 1 = £>%> = £ w </> £p«n = 1 = E^Vi. (23) 

i p q m n 

Now we make use of ER and, for comparison, we use the same basis 
{\a)s>\b) E *} 

\Xim)s\<t>in)E = Y. C aT»S>\b) E ,, \ Vp ) S \%) E = £ D Z\o) S> |&> E ■ (24) 
a, b a,b 

Constraints: 

£ C^Cab H * = ^mm'Snn', E D ab D lb * = <V<W- ( 25 ) 

a, b a,b 



PEi = 22^in\(f>in)E((f>in\, 
n 

(22) 

= 0, while tr E p = £ p Kp\<P P )s{<P P 
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Qp = p-tr E p®p E = E { E A lMn ClrCX"*-E W q D%Dir v }\a)s,{a'\®\b)E>{b' 

a,a',b,b' i,m,n p,q 

(26) 

Hence 

tr E , Qp = 0^ A aa , = E XiPim-KinC'rCTr - E ^DllD p j; = 0, Va, a'. 

i,m,n,b P,<l,b 

(27) 

Again, for a = a': 

E A - = 0, (28) 

a 

as being equivalent with trQp = 0, see eq. (26). 

Validity of eq. (9) assumes validity of eq. (20) for pure and of eq. (27) for 
mixed states. Both eq. (20) and eq. (27) represent the sets of simultaneously 
satisfied equations. We do not claim non-existence of the particular solutions 
to eq. (20) and/or to eq. (27), e.g. for the finite-dimensional systems. 
Nevertheless, we want to emphasize that the number of states they refer to 
is apparently negligible compared to the number of states for which this is 
not the case. For instance, already for the fixed a and a', a small change e.g. 
in ks (while bearing eq. (23) in mind) undermines equality in eq. (27). 

Quantum dynamics is continuous in time. Provided eq. (8) is fulfilled, 
validity of eq. (9) might refer only to a special set of the time instants. So 
we conclude: for the most part of the open S"-system's dynamics, eq. (9) is 
not fulfilled. By exchanging the roles of eq.(8) and eq.(9) above, we obtain 
the reverse conclusion, which completes the proof. 

Q.E.D. 

Lemma 1 establishes: as long as eq. (4) (i.e. eq. (8)) is valid for every 
instant in time, the analogous equality 

ps>(t) = tr E ,p(t) = tr E >Vp(t), (29) 

cannot be fulfilled for the most part of the open ^"-system's dynamics, and 
vice versa. Then, as emphasized above, the projection Qp (Q'p) brings some 
information about the open system S' (S)-in contradiction with the basic 
idea of the Nakajima-Zwanzig projection method. 

The following extends Lemma 1: 
Lemma 2. The two structure-adapted projectors V and V do not mutually 
commute, for every definition (i)-(iii). 

Proof: Let us first consider the projection (i). V'(Vp(t)) = ps'(t) <8> p E >, 
while V(V'p(t)) = p' s {t) ® p' E . Due to eq. (7), p s ,(t) ® p E > ^ p' s (t) <g> pl E ,\/t. 
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The other two projections are mutually similar and assume arbitrary- 
set of the orthogonal operators (states or projectors) for the environment. 
So we consider the case (iii). The projections read: pp{t) = Vp(t) = 
Ei(tr E PEip(t)) <g> P Ei and p' P (t) = V'p(t) = E a (tr E Tl E , a p(t)) <g> U E , a . Then 
V'p P {t) = J2 a (tr E Tl E , a pp(t)) ® Tl E , a = E a crs'a(t) ® H E , a and Vp' P (t) = 
Ei(tr E P Ei pp(t)) <g> P Ei = EiV'si ® Pex- Due to eq. (7): a' Ei (t) <g> P Ei = 
T. m ^m( t )p i s' m ( t )®P i E'm( t )- Hence the commutativity [V,V] = 0, i.e. Vp' P (t) = 
V'pp(t) for every p(t), implies: [p E ' m (t), p E , m ,(t)} = 0, Vi, i', m, m', for ar&z- 
trary choice of the projectors P E i and IT^yq and for every instant in time 
t. This is clearly non-satisfiable. Analogous proof can be given for the case 
(ii)-which completes the proof. 

Q.E.D. 

Thus the Nakajima-Zwanzig projection method faces a limitation. While 
it can be separately performed for any structure (either V or V'), it cannot be 
simultaneously used for a pair of structures-t/ie Nakajima-Zwanzig projection 
method is structure dependent. Once performed, the projection does not in 
general allow for drawing complete information about an alternative struc- 
ture of the composite system-projecting is non-invertible ("irreversible"). 

Our finding refers to all projection-based methods. In formal terms: 
Lemma 1, cf. eq. (29), implies that dVp(t)/dt allows tracing out over 
only one structure of the composite system. If that structure is S + E, 
then tr E 'dVp(t)/dt ^ dp s <{t)/dt [as long as ps*(t) = tr E 'p(t)]. That is, Eq. 
(4) for the S' + E 1 structure is not fulfilled, and therefore cannot provide a 
projection-based master equation for the S' system. This can be seen also 
from the following argument. Tracing out over E is dependent of, but not 
equal to, the tracing out over E'. This dependence follows from the fact that 
the S and E degrees of freedom are intertwined with the S' and E' degrees 
of freedom. Intuitively: "tr E " (e.g. integrating over the £"s degrees of free- 
dom) partly encompasses both the S' and the E' degrees of freedom. On the 
other hand, Lemma 2 excludes simultaneous projecting, i.e. simultaneous 
master equations for the two structures. E.g., dVp(t)/dt = dps{t)/dt <S> p E 
is in conflict with dV'p(t)/dt = dp s >(t)/dt ® p E >: due to eq. (7), only one of 
them can be correct. 

Collecting the told above: Lemma 1 establishes that projection adapted to 
one structure cannot be used for deriving master equation regarding another 
structure, while Lemma 2 emphasizes that simultaneous projecting for the 
two structures is not allowed. Hence, in order to simultaneously describe 
dynamics of the two open systems, S and S', one should avoid projecting of 
the composite system's state. 

Our findings take us back to the beginning, i.e. to Eq. (3), which does 
not have any limitation. It seems that there is not a universal shortcut in 
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deriving master equations regarding the alternative structures of a composite 
system. 

3. Discussion 

We are interested in the variables transformations that simultaneously 
redefine both the open system and its environment. The transformations 
include regrouping of the constituent particles-e.g. in "entanglement swap- 
ping" [13, 21], which is illustrated by eq. (5)-or the more general transfor- 
mations as illustrated by Eq. (6). Such transformations are examples of the 
more general linear canonical transformations performed on the total system 
"system+environment" [19, 20]. The transformations that target only the 
open system without altering the environment can be found in [16, 28]. 

The projections (i)-(iii) in Section 2 are the very special states. Having 
QCR in mind, it is straightforward to show that Lemma 1 of Section 2 equally 
refers to the (ii) and (iii) kinds of projections. 

The time-convolutionless (TCL) method is a particularly useful variant 
of the Nakajima-Zwanzig projection method that employs the projection (i), 
Section 2. The TCL master equation is a time-local differential equation [3] 

^^ = IC(t)Vp(t)+l(t)Qp(t ). (30) 

In Eq. (30) appear the time-local generator K,(t) and the inhomogeneity Z(£), 
while p(t Q ) is the initial state of the total system. 

For the often-used factorizing initial state, p(t a ) = ps(t ) ® p E (t ), the 
projection (i) gives: 

Vp(t ) = p{t ) <£> Qp(t Q ) = 0. (31) 

Eq. (31) directly gives rise to disappearance of the inhomogeneous term 
from Eq. (30) for the S + E structure, and, due to Eq. (7), also for every 
other structure (compare to eq. (9)): 



Vp S (t )®p E (t ) = p S (t )®p E (t ) = p(t ) = V^\iPs>i{to)®PE'i{to) = ^2 XiPS'i(to)<S>PE'i(to) 

i i 

(32) 

However, according to Lemma 1, tracing out the E' environment from 
Vp(t) does not make physical sense. This finding refers even to the mutually 
(macroscopically) "similar" systems S and S'-e.g. if the S' system is obtained 
by joining a single quantum particle (from the environment E) to a multi- 
particle S system (see eq. (5) for a tripartite case). The limitation of the 
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projection methods is a consequence of the quantum correlations relativity 
such as e.g. eq. (7), which is a kinematical result. So, there is not a priori a 
hope that the small changes in the system-environment split would effect in 
small changes in dynamics. 

Encouraged by Eq. (32), one may wonder if there are some sufficiently 
general conditions that can help to circumvent the limitation of the projection 
methods. To the best of our knowledge, the answer should be "no". To 
this end [while bearing in mind exceptions to Lemma 1] the only remaining 
option we can detect is that equality eq. (32) applies to every instant in time. 
However this can be strictly valid only for mutually non-interacting systems. 

Thus we are forced to conclude that the task of simultaneous description 
of the different structures reduces to Eq. (3), yet with the constraint imposed 
by quantum correlations relativity [26]. 

Consider the simplest case for the S+E structure: the tensor product ini- 
tial state, the environment E is harmonic bath of non-interacting oscillators 
weakly interacting with the open system S, applicability of the Born-Markov 
and of the rotating-wave approximation. Then from eq. (7) we can directly 
draw the following conclusions regarding the alternative structure S' + E'\ 
(a) due to the presence of the initial correlations in eq. (7), the S' system's 
dynamics (described by Eq. (3)) is non-Markovian, and also possibly non- 
completely positive [3,4,29-33]; (b) the new environment E' is in general not 
in thermal equilibrium-in general it is in non- stationary state; (c) in general, 
the environment E' consists of mutually interacting particles; (d) one can 
expect non-zero inhomogeneity, X(t), in eq. (30). In addition to this, both 
the strength of interaction and validity of the "rotating wave approximation" 
can be at stake for the alternative S' + E' structure. Thus, in general, inves- 
tigating the alternative open systems dynamics is a formidable task, but see 
[17, 25]. 

Our findings alert: In theoretical considerations, the use of the variables 
transformations in derivation of master equations and interpretations of the 
results should be performed with care. On the other hand, in laboratory 
practice (e.g. in quantum state or process tomography), one should be skiled 
enough to respect locality that is defined by the composite-system's structure 
of interest [19, 20,34]. In regard to the many-particles open systems, which 
in reality dynamically exchange particles with their environments, this can 
pose a non-trivial challenge for experimenter. 

4. Conclusion 

Relativity, i.e. structure dependence, of quantum correlations limits ap- 
plication of the Nakajima-Zwanzig and the related projection methods for 
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certain tasks in open quantum systems theory. Typically, the projection- 
methods-provided information about a subsystem of a composite system in 
an instant of time is insufficient to acquire information about another sub- 
system of the same composite system in the same instant of time. The 
possibility that the two subsystems are macroscopically "almost" the same 
(e.g. in the number of particles) does not in general help. This limitation of 
the projection methods exhibits that the "shortcuts" for describing the open 
systems dynamics may be non-reliable and delicate. 
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